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Abstract. Accurate measurements of the Cosmic Microwave Background (CMB) anisotropy call for high precision and relia- 
bility of the in-flight calibration. For extended surveys the CMB dipole provides an excellent calibration source at frequencies 
lower than ~ 200 GHz; however poorly known foreground emissions, such as diffuse galactic components, complicate the 
signal and introduce a systematic error in the calibration. We show that introducing a weight function that takes into account 
the uncertainty in the a priori knowledge of the sky, allows us to substantially improve the calibration accuracy with respect to 
methods involving galactic latitude cuts. This new method is tested for Planck-LFI radiometers at 30 and 100 GHz. On short 
time scales (less than 1 day) the absolute calibration of each channel can be recovered with an overall 1-2% accuracy. We also 
consider the effect of CMB anisotropy itself on the calibration, and find that knowledge of the CMB pattern on large scales is 
needed to keep the short- time scale calibration accuracy within 1%. 
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1. Introduction 

Accurate determination of the CMB angular power spectrum 
can set powerful constraints on cosmological models and 
allows the determination of main cosmological parameters 
with high precision. Recently, a number of ground-based 
and balloon-borne experiments (IDe Bernardis et aL~2 000 
Hanany et al. 2000| IHaTverson et al. 2002I IPearson et al. 2002l 
Benoit et al. 2003 1 have obtained remarkable evidence for 
the presence of acoustic peaks in the power spectrum, 
extending the pioneering measurements of COBE-DMR 
(Sm ootet al. 1 992) to sub-degree scales, and opening up the 
path to the next generation of space-based precision exper- 
iments ( Bers aneili et al. 2002t . In fact, it has become clear 
(e.g. IDanese et al. 1996b that only space-based observations 
can provide the unique combination of environmental stability, 
freedom from systematic effects and avoidance of ground and 
atmospheric radiation needed for a high-resolution full sky 
survey. The WMAP satellite by NASA was launched on June 
2001 and the first-year data are now available. The Planck 
satellite is an ESA mission planned to be launched in 2007. 

One of the most significant systematic uncertainties af- 
fecting current CMB anisotropy experiments is the instrument 
calibration, typically performed by observing celestial sources 
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with known intensity, such as planets (usually Jupiter), strong 
radio sources (such as Cas A, Carina Nebula, Taurus A) or the 
CMB dipole. The accuracy is limited by the knowledge of the 
absolute flux (~ 5% for planets, > 8 - 10% for radio sources 
and ~ 0.4% for CMB dipole), and also by the limited sky cov- 
erage in the case of the dipole. The WMAP and Planck surveys 
require great calibration accuracy, proportional to their instru- 
ment sensitivities. 

To evaluate the accuracy actually achievable, all possible 
effects (both instrumental and astrophysical) impacting the cal- 
ibration procedure have to be properly evaluated. In this work 
we focus on astrophysical effects by studying the presence of 
different sky emissions contaminating the main CMB dipole 
calibrator. In particular, we propose a method to evaluate the 
accuracy one can expect for absolute calibration on short time 
scales, by using the available information on the sky microwave 
emission. With this method, that can be applied to any precision 
CMB measurement on large sky areas, it will be possible to 
further improve the calibration accuracy any time more precise 
sky observations become available by simply reapplying the 
method we describe. In a forthcoming paper we will analyse in- 
strumental effects (e.g., 1 // noise, pointing uncertainties, ther- 
mal fluctuations, gain drifts) on the calibration performance. 

In the present work we outline the proposed method and 
test the calibration technique for the Planck-LFI receivers. 

The structure of the paper is as follows: in Sect.|2]we briefly 
introduce the concept of calibration; in Sect. [3] dipole con- 
taminations (foregrounds) considered in the analysis are de- 
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scribed. Sect. 0] shows the method to compute statistical and 
systematic errors in the calibration performance, and describes 
the different calibration techniques tested by means of simula- 
tions. Sects. |5]and|6]are devoted to illustrating simulations for 
Planck-LFI and consequent results. In Sect. 0the impact of 
CMB anisotropics on calibration accuracy is studied. The main 
conclusions and future works are discussed in Sect. 03 



2. Calibration concept 

Calibration is the conversion of the output signal of each de- 
tector channel (e.g. telemetry units for a space experiment) to 
physical units (i.e. antenna temperature 1 ). In general, if the in- 
strument response can be considered linear, then 



T A = Go V + T ( 



offset ; 



(1) 



where V is the receiver voltage output corresponding to the sky 
antenna temperature Ta, T^et is an instrumental offset and Go 
is the calibration factor, constant to first order. In practice, Go 
can vary due to instrumental effects, e.g. amplifiers gain or ther- 
mal instabilities. 

For a differential measurement, as in the case of anisotropy 
experiments, the calibration is determined using pairs of 
sources: 



T 2 -T x = G (V 2 - Vi) 



(2) 



where T 2 — T\ - AT a is the antenna temperature difference of 
sources in the sky (often a well-known bright point source and 
the sky background) and V2 - V\ = AV is the corresponding 
radiometer output. 

In order to measure the value of cosmological parameters 
with great precision, calibration must be performed at ~ 1% 
overall accuracy (a simple estimate of calibration requirement 
is shown in AppendixlAl. 

The error in the determination of Go depends on the a priori 
uncertainty on the calibration sources temperatures, ct ATa , and 
on the intrinsic detector noise, cr AV ( Bersan elii et al. 1997i : 



Go 
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°"at a + (Goo-av) 2 



AT A 



(3) 



in this expression cr Co accounts for both the statistical error - in- 
trinsic detector noise - and the systematic error -the uncertainty 
on the temperatures of the calibration sources-. Eq.@ shows 
that a better calibration is performed using higher AT a (as long 
as the corresponding amplitudes do not exceed the linear range 
of detectors). 

Since Eq.[3]is an estimate of the accuracy on Go using only 
a pair of points in the sky, and since in extended surveys much 
more sky pixels are observed, the best value of the calibration 



1 The antenna temperature (T A ) is proportional to the detector 
power P per bandwidth Av: T A = P/(Mv); it depends on both source 
and receiver properties. The relation between antenna and thermody- 
namic temperature is T A = jr^T for intensities and AT A = AT- 



for temperature fluctuations; x = hv/kT, h is the Planck constant, k is 
the Boltzmann constant and v is the observing frequency. 



constant is determined by fitting the distributions of AT A and 
AV, i.e. by minimizing the one parameter^ 2 function: 



X 



AV k - ATf/g 



(4) 



where the index k refers to the pixel pairs available for cal- 
ibration. In the following we indicate with G the value of g 
that minimizes x 2 (g)> while Go is the true value of the calibra- 
tion factor, i.e. the value we need to recover; results will be 
expressed in terms of G/Gq. 

2.1. The CMB Dipole 

An observer in motion with velocity /? = v/c relative to the 
Planckian CMB field sees a dipole pattern: an angular distribu- 
tion of the temperature given by 



°1 -j8cOS(#) 

= T + r oj 0cos(#) + Oifi 2 ) , 



where Tq is the isotropic CMB temperature and § is the an- 
gle between the direction of observation and the direction 
of jS. The CMB dipole has been accurately measured by the 
COBE-FIRAS instrument ( Fixse n et al. 199 6 1 with amplitude 
AT DIP = T /3 = 3.372 ± 0.014 mK (i.e. v = 371 + 1 kms" 1 ) in 
the direction (l,b) = (264.14°±0.15,48.26°±0.15); the overall 
error is 0.4%. 

The dipole is an ideal calibration source for CMB 
anisotropy experiments covering a large sky area. Its amplitude 
is adequate (not too strong to cause non-linear effects) and it 
allows a continuous calibration with no reduction of the obser- 
vation time since it always enters the antenna's field of view. 

3. Foregrounds emission 

Besides the CMB dipole, different components contribute to 
the radio-microwave brightness of the sky. In this context, we 
are interested in these emissions (foregrounds) as "contamina- 
tions" of the prime calibrator, the CMB dipole. In fact, these 
emissions are far less precisely known than the dipole and so 
they represent a drawback for calibration, since they induce 
large errors in the determination of the G factor. In our anal- 
ysis we only considered diffuse components: synchrotron, free 
free and interstellar dust emission. We did not include addi- 
tional diffuse components, such as emission from spinning dust 
grains (Draine & Lazarian 1998 1, since a spatial and spectral 
full-sky template is not yet available for such components. We 
also did not include point sources (e.g. galactic Hn regions or 
supernova remnants) since they fill only a very small fraction 
of the sky pixels. To prove their low impact we considered as 
an example the effect of Hn regions; results will be discussed 
in Sect. El 

Amplitude and spatial distribution of the considered com- 
ponents have been modeled (from available extended surveys) 
to produce synthetic sky maps at frequencies typical of CMB 
anisotropy experiments. We also estimated uncertainties on the 
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intensity of each sky pixel for every foreground component. 
Maps are represented with the HEALPix 2 pixelization scheme 
HGorski et al. 19991 

In this way we are able to produce three sky maps at each 
frequency of interest: a "calibrator" sky, i.e. the best knowledge 
of the sky we can infer today from available data: 

r al («,5) = 2 r - cd(a ' 5) ' (5) 

i 

an "observed" sky, which represents a deviation from T ca i(a, 6) 
according to the uncertainties on the intensity of the considered 
components: 

T ohs (a, T° b \a, S) = J](Tf + cr T ,)(a, 6) (6) 

i 

and a full-sky map of errors cr T (a, 5) obtained adding in 
quadrature the error for each component: 



o- T (a,8)= ^cT 2 Tj ( a ,6). (7) 

Note again that the nature of cr r is systematic. 

The sum in Eq. [5j[6]and0is extended to the components 
considered in the simulation. The terms on the right side of Eq. 
[5]and0are evaluated as follows. 

The systematic uncertainty (1 <x) on the dipole is estimated 
from the COBE uncertainty on its amplitude and direction. To 
be conservative we considered the worst case for calibration 
(see Sect. 16. 1> : when considering dipole, the "observed" com- 
ponent is a "stretched" dipole with v = 370 km s~ 1 , 1 — 263 .99° 
and b = 48.41°. 

To estimate synchrotron emission we used the Haslam et 
al. (1982) all-sky survey at vq = 408 MHz, properly rescaled 
using a spectral indices map: 

T f iv) = Uvo)(^j ■ (8) 

The spatial distribution of the synchrotron spectral in- 
dex jS ; was obtained using the 408 MHz, 1420 MHz 
( |Reich & Reich, 1986) and 2326 MHz dJonas et al. f9 98 1 
maps. To derive errors we took into account the 10% uncer- 
tainty of the original Haslam map plus an error on the spec- 
tral indexes determination; these errors values lie in the range 
[0.1,0.8]. 

Dust emission was estimated using the 100 //m intensity 
and 100/240 fj.m flux ratio map by Schlegel et al. (1998), gen- 
erated from IRAS and COBE/DIRBE data 3 . They also provide 
( Finkbein er et al. 19 99 1 tools for dust emission extrapolation at 
CMB frequencies, using FIRAS data: we used their best model, 
a 2-components model (n.8 in their paper) with mean dust tem- 
peratures < T\ >= 9.4 K and < >= 16.2 K, and spectral 
indices a\ = 1.67 and = 2.70. The authors also estimate the 
errors on extrapolated maps to be ~ 10%. 

2 http : //www . eso . org/ science/healpi x] 

3 Available at http://space.gsfc.nasa.gov/astro/cobe/ 
dirbe_products . html 



Diffuse free-free emission is poorly known. We only 
considered a free-free component spatially correlated with 
dust, with an intensity equal to 30% of dust at 100 GHz 
JDe Zotti et al. 19991 . and rescaled in frequency using a power 
low relation: 

/ y \<*ff 

r(v) = r(ioo)(— ) , (9) 

where v is in GHz and Off = —2.1, the typical Bremmstrahlung 
spectral index. The estimated free-free emission error is 10% 
on dust uncertainties plus 10% on the correlation factor. 

Numerical results for the systematic uncertainty on G due 
to the presence of poorly known foregrounds, highly depend on 
these estimates; anyway, the ability of the weight function tech- 
nique to reduce significantly this uncertainty was successfully 
tested with other scenarios for the "errors" maps, so that the 
proposed method could also be applied in more favourable sit- 
uation, e.g. taking information from the recent WMAP mission 
maps of the microwave emission. 

4. Determination of statistical and systematic 
errors on calibration performance 

To evaluate both the statistical and systematic errors on the re- 
covered gain factor, each complete simulation is split into 2 
steps: 

1 . We set to zero the <tt(q:, 6) terms, so that the same sky map 
is considered: T ohs (a, 6) = T cal (a, 6); this is an ideal situ- 
ation where we assume perfect knowledge of the observed 
signal. This first simulation allows us to evaluate the statis- 
tical component of the error due to instrumental noise. 

2. A sky map T cal {a, S) different from T obs (a, 5) is considered 
(Eq.[5]and|6). In this case both the uncertainties (astrophys- 
ical and instrumental noise) affect the simulation: the dif- 
ference between the result of this second simulation and 
the first one gives the systematic component of the error. 

Of course in both simulations the same noise realization has to 
be used. 

First of all we analyse the ideal case where only the CMB 
dipole is considered; this allows us to evaluate the systematic 
error due to CMB dipole uncertainties. In this simple case, the 
Eq. and are simplified as follows: T ohs = T di P , T cal = 
T d ' p , , . (see SectFJ. 

stretched v * 

Secondly, to consider a more realistic situation, we intro- 
duce galactic emissions: we analyse the impact on G of fore- 
grounds uncertainties only, i.e. in these simulations the same 
dipole is considered both in the observed and in the calibrator 
skies, to fully analyse one type of problem at each time. 

A way to handle the problem of large errors due to fore- 
ground emissions is to cut the galactic plane in the data anal- 
ysis; this was for example the solution adopted by the COBE 
team in the analysis of DMR data ( B ennett et al. 19921 . This 
technique is limited by the presence of emission at high galac- 
tic latitudes, and is not preferable since it forces one not to con- 
sider a portion of data. 
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In order to make maximum use of the whole data set, we 
looked for a technique able to properly weight the signal in- 
tensity T of each point of the sky depending on its uncertainty 
cr T . This can be done introducing a suitable weight function. 
Several functions were considered; based on our analysis we 
concentrated on the family of simple functions 



W(o- T ) 



where 



1 



(To - inf (cr r .) . 



(10) 



(11) 



For every value of a this is a decreasing function, with 
max(W(cr r )) = W(o- ) = 1 (Fig.[T}. 



1.0 

0.8 

^ 0.6 

b 

3, 0A 

0.2 h 



10 



a/cr n 



Fig. 1. Weight function for two fixed values of the a parameter: 
0.75 (solid line) and 1.5 (dashed line). 



When considering pixel-pairs, we use W{ct A t) where 
o- AT = ^cr 2 Ti +crl 2 , 
and 

o-q = inf(cr Ar .) . 
i 

Eq. © is then modified as follows: 



Xw 



(g) = J] W (^T k ) 



AV k -AT^/ g 



°~ 'AV t 



(12) 



(13) 



(14) 



In this way the contribution of pixel-pairs with a poorly 
known temperature is reduced (at a rate controlled by a), lead- 
ing to a great improvement in the calibration accuracy. The 
value of the a parameter to be used is chosen in order to mini- 
mize the overall uncertainty on calibration. 

For a given experiment and scanning strategy, one can study 
the trend of both the statistical and systematic errors, cr stat {a) 
and o- syst (a) respectively, computing the standard deviations of 
calculated errors over the whole mission lifetime as a function 
of a. One expects the weight technique to act in opposite ways 



on the two types of errors: while the systematic error is a de- 
creasing function of a, the statistical error is an increasing one. 
To optimize this procedure one has to search for the minimum 
of the function 

f{a) = [o 2 m {a) + a 2 s J < a)], (15) 
yielding the optimum parameter a. 

5. Calibration for Planck-LFI receivers 

Planck will observe CMB anisotropy and polarization with 
an unprecedent combination of sensitivity, sky coverage, fre- 
quency range and angular resolution. It consists of two instru- 
ments sharing the focal plane of an off-axis aplanatic 1.5 me- 
ter aperture telescope. The Low Frequency Instrument (LFI) 
covers the range 30-100 GHz with four observational chan- 
nels centered at 30, 44, 70 and 100 GHz. The High Frequency 
Instrument (HFI) will observe the sky in six channels between 
100 and 857 GHz with bolometers detectors. The wide fre- 
quency coverage will allow the separation of CMB anisotropies 
from non cosmological signals, thanks to their different spec- 
tral behaviour. 

We applied the absolute calibration procedures described 
in the previous section to the case of Planck-LFI radiome- 
ters. A first analysis of calibration feasibility was performed 
in the context of the COBRAS/SAMBA 4 mission study 
( Bersanel li et al. 1997i . The HFI consortium has analysed the 
problem as well ( Piat et al. 20 02 1, focusing on short time scale 
relative calibration, and long time scale (> 1 month) absolute 
calibration based on the CMB and Earth orbital motion dipoles. 

During operations calibration must be performed with high 
accuracy and as frequently as possible, to control possible 
drifts on instrumental gain and offsets. The monitoring of re- 
ceivers response requires a relative calibration on short time 
scales, while, in principle, absolute calibration can be only per- 
formed at the end of the operations and recovered through the 
whole mission lifetime thanks to the relative calibration. For 
Planck-LFI data, a relative calibration will be performed be- 
tween the 1-hour circles thanks to the observational scanning 
strategy. Regarding the absolute calibration, for LFI receivers 
the main calibration source is the CMB dipole. Only on long 
time scales ( i 3 months) an even better absolute calibration 
can be performed using the dipole modulation due to the satel- 
lite orbital motion around the Sun (which is known with high 
precision), as shown in Piat et al. (2002). On the other hand, 
in the data analysis process it will be of high interest to mon- 
itor the absolute calibration on short time scales (in particular 
for inter-frequency quick-look analysis comparisons) for which 
the dipole modulation cannot be used. The method described in 
this paper allows us to optimise the absolute calibration strat- 
egy of a full-sky mission at short time scales. In addition, this 
concept can be applied to balloon or ground-based experiments 
with limited sky coverage, for which the dipole modulation is 
not an effective calibrator. 



COBRAS/SAMBA was the previous name of the Planck satellite. 
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To simulate LFI receiver calibration, first we created sim- 
ulated data streams 5 considering the Planck baseline scanning 
strategy. According to this strategy, the satellite orbits around 
the L2 Lagrangian point of the Earth-Sun system; the space- 
craft spins at 1 r.p.m. around its spin axis which is kept along 
the ecliptic plane and is repointed by 2.5' every hour to keep the 
anti-solar direction; the telescope field of view is at an angle of 
85° from the spin-axis direction. 

Following the guidelines described in Sect. [3] we gener- 
ate sky maps at Planck-LFI frequencies (30-100 GHz range); 
these maps are then converted into antenna temperature, and 
convolved with a gaussian beam with the Full Width Half 
Maximum (FWHM) of LFI beams (i.e. 33' at 30 GHz and 10' 
at 100 GHz) by the pipeline simulator. The code includes the 
main properties of the Planck payload (e.g. the boresight angle, 
the scanning strategy) and of the considered receiver (e.g. the 
beam location on the focal plane, its FWHM, noise properties). 

The scanning strategy is such that every detector observes 
the same "ring" in the sky for an hour before repointing; in 
our work we average the 60 1 -minute observations of the same 
sky -ring, considering 1 -hour data streams; we refer to these ob- 
servations as "circles". 

The FWHM of the antenna divides every circle ( j) in pixels 
(/). The radiometer output is 



nobs 



+ ri; 



Vij = 



(16) 



where r obs is the observed sky, n is a noise term and Go is 
the true value of the calibration constant. In the noise term we 
only considered white noise, i.e. a random Gaussian distribu- 
tion with rms 



5T n 



nobs 



VAv • ' 



(17) 



where T sys is the system temperature of the receiver, Av is the 
bandwidth (Av/v ~ 20%) and r is the considered integration 
time; the noise amplitude 6T rms is quite constant during the 
mission. In this case, the instrumental uncertainty on output 
differences (see Eq.[3} is: 



O-AV 



4 



5T\ + 5T\ 



Go 



(18) 



where 5T\ and 5T2 are the rms values in the two observed 
points of the sky. 

By averaging 60 rings in one circle we assume that the 
noise of each ring is uncorrelated with the others, which is a 
reasonable assumption for the low 1 // noise expected from the 
instrument (M aino et al. 1999). 

Secondly, we need to determine the pixel-pairs for calibra- 
tion. If the number of pixels in a circle is N, then AT — 1 is 
the number of independent pixel-pairs; we must choose them 
among the N(N - l)/2 possible pairs. Different criteria for 

5 The TOD (Time Ordered Data) used in this work was generated 
by the Planck pipeline simulator of the Level-S of the DPC. Sky simu- 
lations, pointings and other data are available for Planck collaboration 
at http : //planck . mpa-garching . mpg . de/SimData . 



choosing the N — 1 independent pixel pairs were considered: 
the results obtained showed no significant changes. The sim- 
plest approach was then selected: pairs are formed with oppo- 
site pixels in a circle (Fig. [2}: 



1. !6 [l,N/2 + 1] - 

2. is [N/2 + 2,N] 



i' =i + N/2-l 
> V =i-N/2-\ 



In this way every pixel is considered twice; finally the pair with 
the lowest signal difference is not considered, thus leaving ex- 
actly N — 1 pairs. 




Fig. 2. Chosen pixel pairs. 

The k parameter in the^f 2 function (Eq.@} becomes in this 
case k = (m, j), where m = (i, i') refers to pixel-pairs previously 
defined and j refers to the considered circle. 

In our simulations the G factor is computed over two differ- 
ent time-scales: 1 hour and 1 day. Both are "natural" time scales 
for Planck: the first defines the scan "circles" while the second 
is the time scale of Earth-satellite communications. Thus, the 
sum in Eq.@]is extended to the N - I independent pixel pairs 
of one circle in the first case (fixed j value), while is extended 
to the (N— 1) ■ 24 circles covered in one day in the second case. 

6. Simulation results 

We present the analysis carried out at 30 and 100 GHz. The 
first is the channel where foreground emission is the strongest 
while the second one has the highest noise level per single ra- 
diometer (for 30 GHz radiometers T sys ^ 10K, for 100 GHz 
T svs ^ 45K): they represent the worst case for systematic and 
statistical error behaviour respectively. The other LFI frequen- 
cies are in intermediate situations. 

We applied the 2-step procedures described in Sect. 0] as 
already noted, the same noise realization has to be used for both 
steps. In the case of LFI radiometers, since the noise rms value 
has a slight dependence on the observed sky temperature (see 
Eq. I17> . considering the same realization for both simulation 
introduces an error of a few %c on signal terms CEq. I1 6b . 

Similar results hold for all the radiometers at the same fre- 
quency. At 30 GHz we have considered the radiometer LFI-28 



6 



B.Cappellini et al.: Optimized in-flight absolute calibration for CMB surveys 



with the beam position (6b = 5.26°, 4>b = 52.9°) in the foca 
plane and angular resolution of 33' that divides every observe* 
circle in 1950 pixels 5 ; at 100 GHz we used the radiometer LFI 
01 with the beam position {6 B = 2.93°, (f> B = 0.0°) and th 
angular resolution of 10' that divides every observed circle ii 
6498 pixels 6 . 

6.1. CMB Dipole 

As described in Sect.|4l we start by considering the CMB dipol 
only. This allows us to evaluate the systematic error due ti 
CMB dipole uncertainties, and also provides a consistency tes 
of our code. 

Results on the 1-hour time scale are shown in Fig.[3](statis 
tical errors) and Fig. @] (systematic errors); on such short tim 
scales the statistical and the systematic errors are of the sam 
order, i 1- 2%. The typical trend of the statistical error is th^ 

effect of different AT^jp values (<x Co oc AT^ 1 , cfr. Eq. [3} ob- 
served in the scan circles over the different periods of the mis- 
sion time, due to the changing geometry of the field of view 
with respect to the dipole direction. The correlation is clear 
looking at Fig. [6] where the largest temperature difference of 
every circle is plotted. The shape of the systematic error de- 
pends on the choice of v, I and b in the stretched dipole. In Fig. 
@](lower panel) the effect on G when changing these parameters 
one at a time is shown. To be conservative we chose a combina- 
tion (see Sect. of these parameters that maximize the dipole 
systematic error (Fig.|4] upper panel). 








FOB 



1.00 

0.96 
0.92 



















100 200 

Time (days) 



300 



Fig. 3. Statistical error on the G value calculated every hour at 
30 GHz (upper panel) and 100 GHz (lower panel). 

As expected, the amplitude of the statistical and systematic 
errors are consistent with the estimate of detector sensitivities 
on the considered time scale (~ 0.2 mk Hz ^ 2 at 30 GHz and 
~ 0.6 mk Hz~ 1/2 at 100 GHz) and with the CO BE -FIRAS un- 
certainties on dipole measurements. 
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0.02 
0.01 

0.00 
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-0.02 
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100 200 
Time (days) 



300 



Fig. 4. Systematic error due to CMB dipole uncertainty (from 
the COBE measure) on the G value calculated every hour, at 30 
GHz and 100 GHZ (upper panel, plots almost perfectly over- 
lap). In the lower panel, "partial" systematic errors are shown 
(see text for more details); solid line: v = 370 Kms -1 , dotted 
line: v = 372 Kms -1 , dashed line: I = 264.29°, dotted-dashed 
line: / = 263.99°, double dotted-dashed line: b = 48.11°, long 
dashed line: b = 48.41°. 



If calibration is performed every 24 hours, the statistical 
error (Fig. [5} shows the same trend over time but with much 
lower amplitude (cfr. Eq. 1171 : the systematic error is time in- 
dependent and indeed results over different time scales per- 
fectly overlap (results with 24 hours as time-scale are thus not 
shown). These time-trends of statistical and systematic errors 
always apply for a given calibration technique if the only dif- 
ference is the time scale chosen to recover the gain factor. 
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1.000 

0.995 
0.990 

1.010 
1.005 
1.000 

0.995 
0.990 







100 200 

Time (days) 
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Fig. 5. Statistical error on the G value calculated every 24 hours 
at 30 GHz (upper panel) and 100 GHz (lower panel). 



6 Note that the dimension of a pixel is equal to FWHM/3. 
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Fig. 6. Higher AT available for every hour-calibration at 3l 
GHz (solid line) and 100 GHz (dashed line). 



6.2. Foregrounds 

A more realistic situation is obtained introducing galactic emis- 
sions. Histograms for the "error"/"observed" maps of fore- 
grounds (Sect. |5} at the considered frequencies are plotted in 
Fig- El ln b°th cases, for less than 0.4% of the sky pixels this 
value is greater than 1, while isolated pixels have very high 
"error"/"observed "values. 




0.1 1.0 10.0 100.0 0.1 1.0 10.0 100 

error/ signal error/ signal 



Fig. 7. Histograms for the "error"/"observed" maps of fore- 
grounds at 30 (left panel) and 100 (right panel) GHz. 

Fig.|8]shows the systematic error induced by uncertainties 
on foregrounds components. At both frequencies the errors are 
much larger than required (see Appendix lAl and, as expected, 
the situation is far more problematic at 30 GHz. 

The presence of spikes in Fig.|S]means that very large errors 
are concentrated in limited sets of scan circles. In the following 
we give an explanation of the presence of such spikes. 

As already mentioned in Sect. [5] these simulations only 
consider diffuse emission components. Anyway we investi- 
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Fig. 8. Systematic error due to galaxy uncertainty on the G 
value calculated every hour, at 30 GHz (upper panel) and 100 
GHz (lower panel). 



gated the impact of point sources considering galactic Hn 
regions; we used a recent compilation of 1442 sources 
(IPaladini et al. 2003> : they affect less than 0.5% of the sky 
( |Burigana & Pala dini 2003 1, contaminating ~ 7% of scan cir- 
cles in very few pixels. To be conservative, their emission was 
only included in the T obs (a, 6) map: we thus calibrated assum- 
ing no information at all on their presence. Simulations show 
there is no impact on the statistical error, while the systematic 
error is increased (only in the contaminated scan circles): in 
practice other spikes similar to those present in Fig.[8]appear in 
the G-plot. Furthermore we run the calibration procedure ex- 
cluding in the x 2 function (Eq. |4j the pixels with known Hn 
regions; again the statistical error is not affected, but we found 
that most of the spikes in Fig. [S] disappear. This result can be 
explained in a simple way: spikes in the G plots are due to the 
largest values in our estimated error map, which correspond 
-not a surprise- to point source pixels, and to sources not com- 
pletely removed in the maps of the the dust component. 

These tests show that point sources do not interfere signifi- 
cantly with the proposed calibration technique, as long as they 
Ocan be removed with the aid of source catalogues. 

On the other hand strong point sources, such as brightest 
Hn regions and planets, can be very useful as a supplementary 
way to monitor the stability of the gain factor (relative calibra- 
tion), since they are clearly recognizable in the received signal 
given their high intensity. Indeed, relative calibration only re- 
quires stable bright sources, and not a precise knowledge of 
their absolute intensity. When possible, relative calibration on 
these sources can be performed over different time scales, such 
as: 

- less than 1 hour, since they are observed by a given detector 
every minute with a high signal to noise ratio; 

- six months, i.e. the time needed for a given detector to come 
back to the same portion of the sky. 
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6.2.1 . Galactic cuts calibration technique 

As a further step, we simulated the calibration technique with 
a +20° galactic cut; higher galactic cuts do not permit calibra- 
tion over long periods during the mission, due to an excessive 
decrease in the number of avilable pixel-pairs. 

We only show results with a 1-hour integration time, in 
Figs.|5]and[Tni 
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Fig. 9. Statistical error on the G value calculated every hour, 
with a galactic cut of 20°, at 30 GHz (upper panel) and 100 
GHz (lower panel). 
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Fig. 10. Systematic error on the G value calculated every hour, 
with a galactic cut of 20°, at 30 GHz (upper panel) and 100 
GHz (lower panel). The spike in the results is the effect of 
Magellanic Clouds emission at high galactic latitudes. 



As expected, the reduction in the data-set increases the sta- 
tistical error (Fig. [9] compare to Fig.|3J; on the other hand the 



systematic error is highly reduced as can be seen comparing 
Fig-Elk Fig- IH Anyway this technique is limited by the pres- 
ence of residual emission at high galactic latitudes; in Fig. 1101 
we can see the effect of the emission of Magellanic Clouds that 
are at -30° of galactic latitude. At 30 GHz their signals uncer- 
tainty produces an error of ~15% on G. 

6.2.2. The "weight function" calibration technique 

The "weight function" technique allows us to improve the cali- 
bration accuracy. Results on the optimization of a are shown in 
Fig- [HI The optimum choice of a depends on both frequency 
and integration time. The frequency dependency is obvious, 
since foreground emission and noise levels are different at the 
various frequencies. The integration time dependence comes 
from the fact that the systematic error is independent of time, 
while white noise scales as 1 / V?. The a step in simulations is 
0.125 at both frequencies. On the 1-hour time scale, we find 
that the best result is with a = 1.5 at 30 GHz and a = 0.625 
at 100 GHz; on a 24-hour time scale the best result is with 
a = 2.875 at 30 GHz and a = 0.75 at 100 GHz. 

Results for systematic errors are shown in Figs.ll2landll3l 
Considering an integration time of 1 day, one can achieve (T syst 
< 1% for the whole mission time at 30 GHz (Fig.^] upper 
panel) and o- syst is < 0.3% piratically for 100% of the mission 
time at 100 GHz (Fig.[0] lower panel). These results are in line 
with the requirements indicated in AppendixlAl 

7. CMB Anisotropies 

Finally we considered the impact of the presence of CMB 
intrinsic anisotropies on the calibration accuracy. Anisotropy 
maps were simulated with the SYNFAST routine of the 
HEALPix package, given a Standard Cold Dark Matter CMB 
power spectrum. The spectrum was obtained through the 
CMBFAST code 7 . We used the best calibration parameters 
found in our foregrounds analysis. 

The same anisotropy realizations, properly converted into 
antenna temperature, were used at both frequencies; we then 
obtain the observed sky: 



and the calibration sky: 

y^cal _ y^dip , y^gal 



(19) 



(20) 



This is equivalent to neglecting CMB anisotropies in the cal- 
ibration procedure: we assume no information at all on their 
presence. 

Results on systematic errors are shown in the upper panels 
of Figs.[3]andE] Note that the deviations are quite significant, 
exceeding 5% in some periods of the observations. The largest 
deviations are due to large scale CMB anisotropy structures, 
as one can anticipate based on the Planck scanning strategy. 
In fact, if we consider the information on CMB amplitudes on 
large scales (such as the COBE-DMR ones on 9 > 7°) in the 



http : //physics .nyu . edu/matiasz/CMBFAST/cmbf ast .html 
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Fig. 11. Systematic (dot-dashed line) and statistical (dashed 
line) error behaviour as a function of the parameter a (defined 
in Ea. llOl i. The best a value is the one minimizing /(a) (solid 
line, see Ea. ll5t . The horizontal line is the total error for the 
20° galactic cut technique. Panels on the left correspond to a 1 
hour timescale calibration; those on the right to 24 hours. Upper 
panels: 30 GHz channel; lower panels: 100 GHz channel. 

calibrator sky, and LFI-like angular resolution CMB map in the 
observed sky, so that 



j^obs j^dip _|_ j^gal y 



LFI-FWHM 



y^cal y^dip , y^gal , y^cmb 

1 ~ 1 + 1 + 1 COBE-FWHM 



(21) 
(22) 



only small scale structures continue to impact calibration, giv- 
ing an error "randomly" distributed around zero (Fig. 1141 and 
IT31 lower panels). To maintain the systematic error of calibra- 
tion within ~ 1%, it is therefore necessary to use the informa- 
tion on the actual distribution of cold and hot spots in the CMB 
at large scales (e.g. Z 7°). 



8. Conclusions 

The CMB dipole provides a nearly ideal source for calibration 
of extended maps of CMB anisotropy. For sensitive, full-sky 



0.02 
0.01 

o 

£ o.oo 

< 

-0.01 
-0.02 



e5 



7-1 — i — i — i — i — i — i — i — i — | — i — i — i — i — i — i — i — i — i — | — i — i — i — i — i — r 




ry — «y 






100 200 

Time (days) 



300 



Fig. 12. Systematic error on G calculated every hour, with a = 
1.5 at 30 GHz (upper panel) and with a = 0.625 at 100 GHz 
(lower panel). 
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Fig. 13. Systematic error on G calculated every 24 hours, with 
a = 2.875 at 30 GHz (upper panel) and with a = 0.75 at 100 
GHz (lower panel). 



surveys the modulation of the CMB dipole due to Earth ro- 
tation provides the most precisely known absolute calibration 
signal. However, preliminary absolute calibration on short-time 
scales needs to rely on the observed sky signal. Also, experi- 
ments mapping limited sky regions need to use the microwave 
sky as the calibration source. In this paper we have discussed 
an optimisation strategy using the sky emission as an absolute 
calibration source. In the absence of galactic foreground contri- 
butions, which are currently known with rather poor precision 
at millimeter wavelengths, high sensitivity, high resolution ex- 
periments could achieve calibrations accurate at < 1 % level on 
a time scale of 1 hour for nearly the entire mission. However, 
in practice, the presence of galactic diffuse emission introduces 
important systematic effects due to the uncertainties on their 
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Fig. 14. 30 GHz. Upper panel: systematic error on calibration 
due to CMB anisotropics. Lower panel: systematic error calcu- 
lated taking into account information we have on CMB ampli- 
tudes on large (9 > 1°) scales. 
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Fig. 15. 100 GHz. Upper panel: systematic error on calibra- 
tion due to CMB anisotropics. Lower panel: systematic error 
calculated taking into account information we have on CMB 
amplitudes on large (6 > 1") scales. 



intensity and spatial structure. Cuts in galactic latitude can be 
used to mitigate this problem. However, residual systematic de- 
viations in the recovery of the calibration gain are significant 
(typically at the 10% level) due to the presence of high galactic 
latitude structures. We have shown that introducing a suitable 
weight function based on the estimated uncertainty in the to- 
tal signal per pixel, it is possible to improve significantly the 
calibration accuracy with respect to methods simply involving 
cuts in galactic latitude. We applied this concept to the case of 
Planck-LFI, and produced quantitative estimates of the calibra- 



tion accuracy with dedicated simulations. We have computed 
the optimal parameter a for the weight function, which, for a 
given experiment, depends on the frequency and time scale of 
calibration. 

The results show that LFI with its nominal scanning strat- 
egy can reconstruct every hour the value of the gain parameter 
with 1% accuracy for 98% of the total time of the mission at 
30 GHz and for 99.5% at 100 GHz. If we consider an integra- 
tion time of 1 day we obtain an accuracy < 1% for 100% of 
the time at 30 and 100 GHz. Finally we have shown that the 
impact of the presence of the CMB anisotropies themselves is 
less than 1 - 2% for most of the time using the knowledge of 
their large-scale distribution from COBE-DMR. 

While these results have been obtained for the Planck-LFI 
survey, the basic method can be applied to any precision mea- 
surements of the CMB on large sky areas. The recent measure- 
ments at mm wavelenghts by WMAP help to improve the a 
priori knowledge of the galactic signal. These new results can 
be readily incorporated in the proposed technique, leading to 
a more accurate determination of the optimum values of the a 
parameter; we will consider this in a forthcoming paper. 

Future works include the analysis of the impact of instru- 
mental systematics on the calibration procedure; we also need 
to study their behaviour with respect to the weight-function 
technique proposed in this paper. 



Appendix A: Calibration requirement 

The required precision on calibration depends on the goal pre- 
cision in recovering the angular power spectrum coeficients 
C(, typically a few % for precision experiments; this ensures 
cosmological parameters with comparable accuracy. We derive 
here a relation between a systematic error on the calibration 
-the G factor- and the accuracy on the recovered power spec- 
trum C{, in the simple case where G is constant over the whole 
map. 

An error 6G on the gain factor induces an error on the mea- 
sured AT: 



AT = G AV 
6G 



(> + t) 



AV : 



Ar 



then 



5G 

AT -AT= — Ar. 
G 



(A.l) 



(A.2) 



In the first case the coefficients of the spherical decomposition 
of the Ar function are 



Ar 



(e,<t>)Y (m (e,<p)d£i, 



while in the second one 



r at' 

a (m - J —(0, 



</>)Y (m (e,(f>)dn. 



(A3) 



(A.4) 
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Subtracting the previous two expression, 

r AT' - AT 
Sa {m = I (6,<t>)Y tm (6,<t>)dQ. 



I 



T 

SG AT 

— —{6,4>)YtJ0,WLn 



5G 

= -pra (m . (A.5) 
We recall the relation between Q and af m coefficients: 

1 1 

c t = 277Y X |flft « |2 ; (A - 6) 

m=-t 

differentiating 

dC 2 - 

5C ( = 5a Cm = — — - } \a tm \5a tm 

da (m 2t + 1 £j ( 

2 5G -A , _6G 



^ \a tm \ 2 = 2— Q (A.7) 



m=- 



21+ 1 G 
5C C 5G 

-» 7^ = 2— - (A.8) 
In conclusion, in order to get 6C(/Ce ^ 1 - 2% one requires 
^<1%. (A.9) 
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